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ABSTRACT 

A  harmonic  analysis  of  the  observed  diurnal  variation  of  the 
coefficient  of  eddy  viscosity  is  conducted  at  selected  levels  followed 
by  a  mathematical  analysis  of  the  amplitude  coefficients  in  order  to 
determine  a  function  of  height  and  time  needed  to  represent  this 
variation,.   It  is  found  by  this  approach  that  a  lengthy  series  of  terms,,, 
consisting  of  exponential  and  polynomial  functions  of  height  and 
trigonometric  functions  of  time  and  height,  would  be  needed  to  represent 
the  observed  variation,, 

A  two  layer  empirical  model  is  developed  from  the  observed  data 
which  gives  the  coefficient  of  eddy  viscosity  as  a  function  of  stability 
and  height  in  the  lower  layer ,  and  as  a  function  of  vertical  shear  of  the 
horizontal  wind  and  height  in  the  upper  layer,, 

The  author  wishes  to  express  his  appreciation  for  the  assistance  and 
encouragement  of  Professors  G0  J»  Haltiner  and  A„  B„  Mewborn  of  the 
Uo  So  Naval  Postgraduate  School  in  this  investigation,, 
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1,    Introduction. 

In  an  attempt  to  account  for  the  observed  wind  structure  in  the 
friction  layer,  meteorologists  have  devoted  considerable  effort  to 
determine  the  variation  of  the  coefficient  of  eddy  viscosity  with 
respect  to  time  and  space* 

The  layer  of  frictional  influence  has  been  subdivided  into  the 
laminar,  surface,  and  spiral  layers  each  possessing  certain 
characteristic  properties.   The  coefficient  of  eddy  viscosity,  M      , 
has  been  shown  to  vary  linearly  with  height  in  the  neutral  surface 
layer  and  as  a  power  function  of  height  in  the  non-neutral  surface 
layer.   Using  wind  data  obtained  from  double  theodolite  observations 
and  assuming  a  steady  state,  Kildner  _6 J  solved  the  horizontal 
equation  of  motion  for  M     utilizing  finite  differences  for  the 
derivatives  of  velocity.   He  showed  that  ^M     increased  up  to  a  fairly 
low  elevation  in  the  spiral  layer  (2^0  meters)  and  then  decreased  to 
a  small  residual  value  near  the  gradient  level. 

The  Taylor   7  |  solution  for  the  wind  in  the  spiral  iayeT  assumes 
fA      constant,  along  with  a  number  of  other  broad  assumptions,  and 
yields  a  general  representation  of  the  observed  vertical  wind  variation, 
More  recently,  Buajitti  and  Blackadar  |_  2j  and  Haltiner  I  4 J  have 
sought  the  distribution  of  eddy  viscosity  as  a  function  of  time  and 
height  in  order  to  account  for  the  diurnal  variations  of  the  wind 
structure  in  the  friction  layer.   Haltiner  [  4_|  assumes  the 
coefficient  of  eddy  viscosity  to  have  the  form 

Ki/f   =  g(t)  k(z) 


v/iiere 


g(t)  =  1  +  b  (Sin  Co  t  +  c  Sin  2co  t) 

and  k(z)  takes  on  several  different  exponential  forms o 

It  is  the  purpose  of  this  investigation  to  determine  the  degree 
of  complexity  of  the  function  of  /A     (z,  t)  needed  to  yield  observed 
values  of  yM  ;  and  also  to  determine,  if  possible.,  the  coefficient 
of  eddy  viscosity  as  a  function  of  stability  and  vertical  shear  of  the 
wind. 


2o   The  Coefficient  of  Eddy  Viscosity  as  a  Function  of  Height  and  Time0 

The  horizontal  equations  of  motion  for  turbulent  flov;  where  the  eddy 
stresses  are  assumed  to  be  a  function  of  height  and  time  only  and  the 
molecular  or  viscous  stress  is  considered  negligible  are 

bt  r  ^ 

Here  l^K  an^  '  ZM  are  ^e  components  of  eddy  stressj 
^q  and  *V"q  the  components  of  the  geostrophic  windj  V   ,  the  mean 

wind  velocity  vector  of  which  U  and"V  are  the  horizontal  components; 
p   ,  the  density;  and  -T   ,  the  coriolis  parameter,, 

The  eddy  stresses  are  assumed  to  be  related  to  the  mean  wind  shear 

by  the  expressions 

bz  t)z  U) 

The  geostrophic  departure  method  employs  observations  of  the  mean 
velocity  and  geostrophic  wind  to  evaluate  all  the  terms  of  eq  (1),  so 
as  to  yield  the  distribution  of  ^  ^*    and  ^!  -  z3   0   This  method 
gives  the  distribution  of  stress  and  eddy  viscosity  provided  the  wind 
components  have  a  maximum  at  some  levels  h  and  hy  so  that  at  hx,xr- --V 
and  at  h^,.  --;0  . 

J    oz. 

By  multiplying  eq  (1)  by  ^P   and  integrating  from  hx  to  an 
arbitrary  height  z,  the  distribution  of  J^x  is  obtainedo   Similarly 
the  distribution  of  lzu  is  foundo   Then  eq.(2)  leads  to  two  independent 
estimates  of  the  distribution  of  eddy  viscosity,, 


Assumin 


g  no  discrepancies  in  the  procedure  or  data5  ftx-Mu^O. 


Blackadar  f~lj     used  the  geostrophic  departure  method  on  a 
composite  wind  distribution  obtained  by  averaging  the  winds  at  two 
stations  (Oklahoma  City,  Oklahoma  and  Wichita,  Kansas)  over  24 
selected,  not  necessarily  consecutive,  one-day  periods „   Figures  1 
and  2,  reproduced  from  Blackadar  9s  £lj  data,  show  the  eastward 
and  northward  components  of  the  average  winds  during  the  24  day 
period 0   Although  great  care  was  taken  in  determining  each  of  the 
terms  of  eq  (l),  Blackadar  '  [_lj  pointed  out  that 


"the  independent  estimates  of  J^x    and/^  did  not 
agree  and  negative  values  of  each  were  numerous  and 
frequently  large „   It  is  assumed  that  the  reason  for 
the  poor  performance  is  in  the  inaccuracy  of  the 
geostrophic  windD " 


For  this  reason,  he  devised  an  indirect  method  for  obtaining  the 
values  of  Mq  and-Vq  which  would  minimize  the  sum  of  the  absolute 
magnitudes  of'  /AK-Mu   at  the  four  levels  (950,  900*  8^0  and  800  mb)  „ 
He  forced  MK  to  be  equal  to  Ms*    at  a  single  level  and  at  the  two 
standard  radiosonde  times,*  then  the  geostrophic  wind  components  were 
computedo   Next  utilizing  the  thermal  wind  relationship  and  assuming 
the  temperature  sounding  to  be  accurate,  the  values  of  Mq  and'Va  were 
found  in  the  vertical.   The  variation  of  geostrophic  wind  at  times 
other  than  the  standard  radiosonde  times  was  determined  on  the 
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i?he  eddy 

or  height 

intervals  and  three  ho  The  values  of     and 

so  obtained  were  general;  at  all  times,,   Blackadar   1 

deca  ded 

"oooto  dispose  of  the  differences.,  the  beat 

estimate      as  assumed  to  be  given  by 

Figo  5  shows  Blackadar 3s   I      aputed  values  of     as  a  function  of 
height  and  time0 

It  is  desired  to  find  a       oh     :'z5,t)  which  will  yield  these 
values  of     „   Due  to  the        c  nature  of  the  solution,  figo  3 
must  be  modified  in  order  to  show  the  sane  values  of 
corresponding  heights  at  the  beginning  and  end  of  the  24  hour  periods, 
The  alteration  of -fig.  3  was  minimized  insofar  as  possible,  arid  the 
modified  figo  5  is  shown  as  fig„ 

The  following  approa<     the  problem  was  used,,  Determine  the 
time  series  which  fits  the  observed  values  of  /-:  for  12  times 
(  2  hour  intervals)  at  each  of  the  given  heights  0,  200 s.  400,  600oOOO 
l£00  meter So  An  example  for  the  height  of  200  meters  is 


(4) 


2  m/A  (2,t)  -  Ac  (2)  4  A1  (2)  cos   [  (t  +  Tj  (2)  J 
A?  (2)  cos  [(2i  -   (2) J  +,aook6   cos  £(6t  +  T$(2)] 


Kere  sine  and  cosine  terms  belonging  to  the  same  harmonic  have  been 
combined'.  Aq(2)  is  the  mean  value,  A]  (2)  is  the  amplitude  of  the 
first  harmonic,  T|(2)  is  th  Le  of  the  first  harmonic  in 

degrees  and  t  is  time  in  hours  after  t  ~  00  (1500  CST),  expressed  in 
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decrees,   Now,  having  10  equations  of  the  same  form  as  eq  (4-),  fit  a 
continuous  function  of  z  to  the  10  values  of  A  and  T  for  each  harmonic 
which  makes  a  significant  contribution,,   This  would  yield  the  final 
form  of  the  equation, 

M   (z,t)=  AQ(z)  +  ^ (z)  cos  [t  +  Tx(z)  I 
+  A2(z)  cos  |_2t  +  T2(zv,    +  . ....A£(z)  cos   |6t  +  r^(z)J 

The  harmonic  analysis  at  each  even  200  meter  height  was  accomplished 
yielding  the  values  in  Table  1. 

Fig.  5  shows  values  of  Aq  (the  mean),  Ai  and  A2  versus  height; 
fig,  6  shows  A3,  A4  and  Ak  versus  height.   Fig.  7  shows  T^  and  T2 
versus  height  and  fig.  8  shows  T^  TA_  and  Tk  versus  height. 

A  diagonal  difference  table  was  computed  for  Aq  values  versus 

2 

height.   The  second  differences  of  Aq  &   Aq^ showed  the  most  nearly 

constant  values.  Although  the  second  differences  were  not  very  close 
(ranging  from  4^  to  128),  a  parabola  was  fitted  to  the  values  of  Aq 
by  the  least  squares  method.  As  expected,  the  parabola  did  not  fit  the 
data  to  within  sufficient  accuracy. 

Next  it  was  noted  from  figs.  5  a^d  6  that  a  possible  orderly 
decrease  in  the  ratios  of  A./Aq  existed  at  each  height  a;,  i  increased 
from  1  to  5»   A  plot  of  the  ratios  A-/Aq  versus  height  was  drawn 
(see  fig.  9).   Utilizing  these  ratios  the  equation  for  M  (z,t)  would 
take  on  the  form 

JA   =  A0(z)  £l  +BX  (z)  cos  Tt  ♦  T1(z)1   +'B2(z)  cos  T(2t  4  T2(z)] 

+.-... Bg  (z)  cos  \{6t   +  T6  (Z;j} 
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TABLE  1. 

Results  of  harmonic  analysis  at  10  levels,, 

Z  is  height  above  surface  (m) ,  A~  the  mean  coefficient  of  eddy  viscosity 
(g  cm"  sec"  )y  A^  the  amplitude  of  the  i'tn  harmonic  (g  cnT  sec  )  and 
Ti  is  phase  angle  of  i   harmonic  degrees„  ^  ' 


z 

Ao 

h 

Tl 

A2 

T2 

A5 

T5 

A4 

T4 

A5 

T5 

A6 

0 

204.1 

255  0  4 

26.1 

86.2 

40.2 

4o.O 

58.7 

22.9 

54.9 

7.6 

96 .8 

2.6 

200 

118.5 

158.8 

19.7 

74.8 

54.9 

55o6 

61.2 

24.4 

85.7 

19.4 

125.8 

8.8 

4oo 

75.8 

87.8 

12.9 

54.5 

51.9 

24.8 

69.1 

16.6 

115.5 

18. 5 

156.0 

9.5 

6oo 

87.0 

105.5 

15.7 

50.0 

54.5 

19.4 

69.4 

15.5 

125.8 

25.7 

I65.I 

12.2 

800 

144.5 

191.9 

11.9 

98.7 

27.6 

45.8 

45.7 

18.5 

106.5 

26.0 

152.2 

14.4 

1000 

266 

577 « 6 

7.7 

255.5 

16.0 

158.1 

18.6 

95.0 

19.5 

68.6 

18.7 

51.2 

1200 

461.2 

690.O 

5.5 

471.5 

8.2 

551.1 

8.5 

248.4 

6.5 

196.6 

6.9 

95.7 

l4oo 

758 

II62 

2.2 

806 

7.1 

560 

8.5 

560 

8.0 

258 

8.0 

114.6 

1600 

1160 

1761 

2.1 

1169 

8.6 

756 

11.8 

458 

15.6 

266 

15.2 

105 

1800 

1690 

2518 

5.4 

1620 

18.1 

984 

51.0 

559 

44.0 

249 

45.5 

60 

(1)  All  surface  values  are  the  result  of  extrapolating  Blackadar's  j^lj 
curves  down  to  the  surface. 
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Here  Aq(z)  is  the  function  of  height  needed  to  fit  the  observed  values 
of  the  mean  M  at  each  height,   %(z)  is  the  function  of  height  needed 
to  fit  the  ratio  A^/Aq  at  each  height  and  T^(z)  is  the  function  of 
height  needed  to  fit  the  observed  Tj_  values  at  each  height.   If  the 
curves  of  fig,  9  could  be  fitted  by  a  relatively  simple  function  of 
height  vice  an  estimated  5*h  °f   6th  degree  polynomial s  there  existed  a 
possibility  of  a  series  analysis  and  possible  further  simplification. 
One  noteworthy  aspect  of  this  approach  should  be  observed.   It  appears 
that  a  fairly  simple  function  of  height  (e»g,,  a  parabola)  could  be 
fitted  to  the  ratio  curves  (fig,  9)  with  reasonable  accuracy  if  a 
layerized  approach  is  followed.   The  dashed  lines  on  fig,  9  represent 
the  probable  curves  which  would  fit  a  2  layer  model.   The  division  of 
layers  appears  to  be  in  the  vicinity  of  600  to  700  meters.   This 
conceivably  could  be  of  physical  significance. 

Next  it  was  noted  that  Aq  (fig,  5)  increases  non~linearly  with 
height  above  46 0  meters.   Further,  the  point  z  =  46 0  meters ,  for  which 
A~  =  75g.3m~  sec"  „  appears  to  be  a  point  of  symmetry.   The  axis  was 
translated  to  this  point  and  Aq  versus  height  was  plotted  on  log-log 
paper;  this  plot  coincided  almost  exactly  with  a  straight  line. 
By  the  least  squares  method,,  Aq(z)  =  az  <*  b  was  determined. 


The  equation  for  A-.  was 


2,1 


AQ(z)  -  6,C     z  -  4,6  p1  *  75  (7) 

where  z  is  height  above  the  ground  in  hundreds  of  meters.   Table  2 
shows  a  comparison  of  the  computed  values  of  Aq  and  observed  values 
(from  the  harmonic  analysis). 
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TABLE  2. 
A  (z)(g  cm  sec"  )  computed  versus  A~(z)  observed, 


HEIGHT 

Z  .    x  10~2 
meters 

COMPUTED 
A0(Z) 

OBSERVED 
A0(Z) 

PERCENT 

ERROR 

0 

225 

204„ 

10.2 

2 

118. 

118. 

0.5 

4 

75 

75< 

5.6 

6 

35- 

37 

2.0 

8 

155- 

144. 

6.5 

10 

286. 

266 

7.5 

12 

487. 

461, 

5.7 

14 

759* 

758 

0.2 

16 

1105 

1160 

4.7 

15 

1524 

1690 

9.8 
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Similarly  A,  (z)  and  A?(z)  were  obtained.  A  comparison  of  A-^(z) 
computed  and  observed  is  shown  in  Table  5>  an(i  for  A~(z)  in  Table  4. 
The  equations  developed  by  the  least  squares  method  are  as  follows: 


Ax(z)  =  8.4 


-  4.' 


2.15 


+  81 


(8) 


A0(z)  =  2.5 


-  A.o 


2.4b 


4  55 


(9) 


When  locating  the  approximate  point  of  symmetry  of  the  Ap  versus 
height  curve  (fig.  5) »  it  was  necessary  to  employ  the  method  of  least 
squares  to  determine  which  of  three  logically  chosen  positions  for  the 
height  axis  yielded  the  least  difference  between  A  (z)  calculated  and 
Ap(z)  observed.   In  spite  of  this  effort,  relatively  large  differences 
were  observed  in  the  range,  surface  to  400  meters  (see  table  4),  v/here 
the  symmetry  was  weakest;  however,  the  fit  appears  to   e  good  above 
oOC  meters. 

A  linear  approximation  was  used  to  fit  T-,(z)  and  T0(z)  (see  fig.  7). 
The  equations  determined  by  the  least  squares  method  are 


T,(z)  =  -  1.20  z  +  21.4- 


T0(z)  =  -  1.81  z  +  59 


(10) 


where  z  =  height  above  surface  in  hundreds  of  peters  and  T-,  (z)  and  T~(z) 
are  in  degrees.   Table  5  shows  a  comparison  of  the  phase  angles  computed 
1  y  eq  (10)  for  the  first  and  second  harmonics  and  the  observed  phase 
angles. 
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TABLE  5. 
A-i(z)(g  :m""  sec"  )  computed  versus  A^(z)  observed 


HEIGHT 
Meters  x   lcP 

COMPUTED 
Ai(Z) 

OBSERVED 
AX(Z) 

PERCENT 

ERROR 

0 

295. 

255. 

25.9 

2 

l4l. 

Hi 

1.7 

4 

81. 

88. 

7.6 

6 

105 

105. 

0.5 

8 

205 

194 

6.8 

10 

408. 

577. 

8.2 

12 

719. 

690 

4.5 

l4 

1142 

1162 

1.7 

16 

1685 

1761 

4.5 

18 

254o 

2518 

7.1 
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TABLE  4. 


-1   -1 


A?(z)(g  3m~  sec   )  computed  versus  A2(z)  observed 


HEIGHT 
^neters  x   *-® 

COMPUTED 
A2(Z) 

OBSERVED 
AX(Z) 

PERCENT 
ERROR 

0 

110. 

86. 

27.6 

2 

48. 

74. 

54.9 

4 

55- 

54. 

55.5 

6 

48. 

50. 

2.6 

8 

110. 

98. 

11.8 

10 

27P 

255. 

2.4 

12 

kSO. 

471. 

4.8 

14 

751. 

806 

6.8 

16 

1157 

II69 

1.0 

18 

16?4 

1620 

5.5 
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TABLE  5. 

T, (z)  computed  versus  T^(z)  observed 

and  To(z)  computed  versus  ?2(z;  observed  for  10  levels. 

T,(z)  and  T0(z)  in  degrees. 


-2 
Z     .           x   10 
neters 

T1(z)obs. 

T]_(z)conp. 

T2'z)obs. 

T2(z)conp. 

0 

26. 

°1 

4o. 

59 

2 

19. 

18. 

5£ 

55. 

4 

13 

16.5 

5  & 

3* 

6 

1^ 

14.1 

5*. 

28. 

8 

3ft. 

1* 

27. 

24. 

10 

7. 

9- 

16. 

20. 

12 

5- 

7 

8. 

17. 

14 

2. 

4. 

7. 

15- 

16 

2. 

2.0 

8. 

10. 

18 

5- 

-0.4 

IS. 

6. 
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in  examination  of  the  remaining  amplitude  coefficients  and  their 
corresponding  phase  angles  (see  figs.  6  and  8),  it  appeared  that  the 
functions  of  height  needed  to  fit  these  values  to  the  desired  degree  of 
accuracy  would  be  relatively  long  and  too  complex  to  warrant  further 
effort  in  this  phase. 

Also  one  of  the  objectives  was  to  see  how  complicated  a  function  of 
height  and  time  it  would  take  to  yield  the  observed  values  of  yM  in  fig.  3' 
The  equation  foryW  (z,t)  thus  far  obtained  is 

J<A  (z,t)  =  6.0   Is  -  4.6  |2J  +72+:G.£f   z  -  4.5  j  2'  %   +  8l] 


cos 


COS 


L  24 


+  (21 


-  1.2 


<0^0>  [2.^1 


lso j'  r 


-  4 


2.5~ 


4  55 


'1 


(11) 


1  2irt   /,0   ,  q  *  ii  "I 


.-? 


where  z  =  height  above  ground  in  meters  x  lO"*"  and  t  is  time  in  hours 
after  t  =  0  (1^00  GST).   Values  of  M  (z,t)  were  computed  using  the 
liar  chant  desk  calculator  for  time  intervals  of  one  hour  and  every  100 
meters  height.   Fig.  10  shows  computed  values  of  /-/plotted  on  a 
height-time  coordinate  system. 

Comparison  of  fig.  10  and  the  desired^  distribution  from  fig.  4, 
shows  fig.  10  to  contain  a  fairly  representative  set  of  isolines  for  H 
from  1000  C3T  to  2000  C37  from  the  surface  to  1800  meters.   Both 
figures  show  the  largest  order  of  magnitude  for  the  coefficient  of 
eddy  viscosity  at  approximately  1500  C3T  at  all  levels.   In  addition, 
both  figures  show  a  minimum  f or /^  in  the  vertical  at  approximately 
400  meters  between  1000  C3T  and  2000  C3T  withyM  decreasing  from  the 
surface  to  400  meters  throughout  the  period  and  with  ft  increasing 
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rapidly  with  height  above  400  meters  from  1200  CST  to  1300  OST.   In 
fig.  10,  however,  the  decrease  of  /4    with  time  at  any  level  above  1000 
meters  is  too  great  in  the  vicinity  of  1900  OST,  as  is  the  increase  of 
yM  with  time  too  large  at  any  level  above  1000  meters  in  the  vicinity 
of  1000  C3T.   Also,  the  vertical  gradient  of  M  above  400  meters  at 
1^00  CST  is  less  than  that  observed  on  fig.  4.  ±he   latter  could  be 
rectified  by  computing  the  remaining  harmonics  which  would  tend  to 
increase  the  values  of  fi  computed  at  1^00  OST  (t  =  0)  so  that  they 
more  closely  approach  corresponding  values  of  observed/^  (fig.  10). 

Between  2000  OST  and  1000  C3T,  surface  to  1800  meters,  a 
comparison  of  fig.  h   and  fig.  10  shows  fig.  10  to  have  a  completely 
unacceptable  representation  of  the  distribution  of jM  .      Fig.  10  shows 
a  column  of  minimum  M  values  extending  vertically  from  the  surface  to 
1800  meters  at  about  2200  CST,  with  no  phase  lag,  and  another  minimum 
column,  extending  from  the  surface,  at  0^00  CST,  to  1800  meters,  where 
the  minimum  occurs  at  about  0700  CST.   These  minima  are  separated  by 
a  relatively  strong  maximum,  surface  to  1800  meters,  occurring  at 
0200  CST.   In  contrast,  the  desired  form,  fig.  4,  shows  a  strong 
single  minimum  extending  from  the  surface,  where  it  is  estimated  to 
occur  at  2^00  CST,  to  1800  meters  where  it  is  observed  at  approximately 
0500  CST. 

The  double  maxima,  double  minima  vice  the  desired  single  maximum, 
single  minimum,  are  attributed  primarily  to  the  fact  that  the  mean  and 
only  two  of  the  six  harmonics  were  utilized  in  deriving  the  values  of 
fA.   for  fig.  10.   Fig.  11  is  introduced  to  lend  credence  to  the  above 
remark.   This  figure  shows  the  resultant  of  the  first  two  harmonics 
versus  time  for  height  equal  to  1200  meters.   The  resultant  M 
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distribution  utilizing  the  mean  and  the  first  two  harmonics  as 
determined  by  the  harmonic  analysis,  not  from  fitted  values ,  tends  to 
show  the  following  characteristicsSan  absolute  maximum  at  1500  C3T,  a 
relative  maximum  at  0250  GST,  and  two  minima,  one  at  0700  03T  and  the 
other  at  2200  C3T0 

Each  of  the  above  maxima  and  minima  are  observed  on  fig.  10  at 
the  corresponding  times  and  the  magnitude  comparison  is  as  follows; 

TIME  ^    from  fig.  10         yCf  from  fig,  11 

1500  C3T  I65O  1620 

2200  C3T  -85  -l4o 

02^0  C3T  215  240 

0700  GST  -124  -140 

Similar  analysis  at  other  heights  would  yield  similar  resultSo 

This  would  tend  to  show  then  that  the  fitted  values  have  not 
distorted  the  picture  of  At  distribution  with  time  and  height  except  to 
a  slight  degree  and  that  the  cause  of  the  deviation  of  the  computed  JA 
values  (fig.  10)  from  21ackadarss  ML  j  observed  values  (altered 
fig«  4)  must  be  the  remaining  four  harmonics .   Further.,  when  using 
eo  (11)  and  including  the  remaining  four  harmonics,  the  values  of 
M    computed  compare  very  favorably  with  the  observed  values,, 

Pertinent  to  the  variation  of  M  with  height,  Estoque  f5  | 
has  recently  approximated  the  distribution  of  M   by  a  linear  decrease 
ofM   above  a  maximum  A4  value  found  at  a  height  of  about  50  meters,, 
He  proposes  a  two  layer  model  for  the  frictional  layer,  the  lower 
layer,  called  the  surface  sub-layer,  extends  to  approximately  50 
meters,  and  an  upper  layer,  called  the  transition  zone,  extends  to 
approximately  2000  meters „   Utilizing  this  fi  distribution  he  has 
numerically  integrated  the  equations  of  motion  and  hect  conduction 
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equation  bo  obtain  the  24  hour  distribution  of  vri.net  and  temperature. 
The  diurnal  wind  distribution  thus  obtained  by  Sstoque  I  5  J  appears 
to  be  in  fairly  good  agreement  with  the  diurnal  wind  characteristics 
observed  on  several  occasions „   However,  for  Dlackadar's  |_1  J  data, 
(see  figo  5)s  it  was  necessary  to  attain  an  average  height  of 
approximately  ll^O  meters  before  Ja  varied  linearly  (but  increased) 
with  height  to  within  a  10  percent  error  of  the  observed  value  of  M  „ 
Further,  the  slopes  of  the  lines  representing  the  linear  variation  of 
j^K   with  height  (for  different  times)  varied  up  to  55°» 
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Vertical  Shear  of  the  Wind. 

Wind  data  for  this  phase  were  obtained  from.  figs.  1  and  2, 
reproduced  from  Blackadar  |_1_|   .   Temperature  data  for  the  24  selected 
days  vere  obtained  fron  photo  copies  of  WBAN-55  f°r  *he  months  of  June, 
July,  August,  and  September  1952  for  Wichita,  Kansas  and  Oklahoma  City, 
Oklahoma,  which  vere  furnished  by  the  national  Weather  Records  Center, 
Ashville,  North  Carolina,, 

All  temperatures  at  the  reported  levels  (surface,  950>  900,  850, 
800  and  750  mb)  at  both  cities  and  at  both  sounding  times  (1500  GOT 
and  0^00  OCT)  were  converted  to  potential  temperature.   The  surface 
potential  temperature  was  obtained  by  the  Arowagram  and  upper  level 
potential  temperature  uas  obtained  from  the  Smithsonian  Meteorological 
Tables. 

For  the  0J00  OCT  sounding  and  at  each  level,  the  average  potential 
temperature  uas  obtained  by  adding  the  24  observations  at  Wichita  to 
the  24  observations  at  Oklahoma  City  and  dividing  by  48.  This  yielded 
the  average  sounding  for  the  two  cities  over  the  24  day  period  at  0^00 
OCT.  A  similar  procedure  was  followed  to  obtain  the  average  height  of 
the  950,  900,  850,  800  and  750  mb  levels.  All  heights  are  in  terms  of 
O.98  geo dynamic  meters.   The  same  procedure  was  followed  at  1500  OCT. 

The  average  heights  and  potential  temperatures  obtained  for  the 
two  times  are  shown  in  table  6.   Next  the  stability  factor 

e  bz 

was  computed  for  the  center  of  the  averaged  layers  using  finite 
differences  for  the  partial  derivative  of  &       .   Stability  factor 
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TABLE  6. 

Twenty  four  day  average  heights  (in  0.98  geo dynamic  meters) 
and  potential  temperatures  (degrees  K)  for  the  reported 
levels  of  concern. 


0^00  GOT 

1500  GGT 

Level 

Z 

e 

Level 

Z 

9 

meters 

degrees  K 

meters 

degrees  K 

Surface 

0 

500  oS 

Surface 

0 

501.8 

950mb 

188 

505.0 

950mb 

201 

502.5 

900mb 

659.1 

506,3 

90Cmb 

672.5 

505.1 

S^Omb 

1155.5 

507  j. 

850mb 

1168.2 

507  oC 

800mb 

1677 

509.3 

SOOmb 

1687.9 

508.  b 

750mb 

2228.2 

5H.5 

750mb 

2257.6 

510.7 
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versus  height  of  center  of  layer  for  the  two  times  is  shown  in  table  7° 

Values  of  the  components  of  the  vertical  wind  shear,  ^ —  and  STT  > 
for  each  100  meter  level  were  then  obtained  from  figs.  1  and  2 
utilizing  finite  differences  for  the  partial  derivatives .   The  shear 
factor,  J  (  tr —  )  "W  - —   //^  „   \/as  computed  next  for  each  100  meter 
level  for  the  tuo  times. 

Graphs  of  M  observed  versus  the  stability  factor  and  of/M observed 
versus  the  shear  factor  uere  then  constructed  in  an  attempt  to  find  a 
correlation  between  these  factors  and  M   .   One  apparent  correlation  was 
that  in  the  lower  layers  (below  600  meters)  where  pronounced  stability 
variations  occurred,  strong  stability  was  associated  with  a  small  M    , 
and  vice  versa.   Above  about  600  meters  the  contrast  in  stability  factor: 
was  considered  negligible.  Moreover,  there  was  insufficient  difference 
between  the  M  values  at  each  level  for  these  two  times  (0900  C5T  and 
2100  G3T) ,  nor  between  the  shear  factors  in  the  upper  levels  (above 
600  meters),  to  detect  any  correlation  between  wind  shear  and  the 
values  of/1  observed.   In  the  lower  layers,  when  the  shear  appeared  to 
be  large,  the  value  of  M  was  also  large.   In  order  to  see  if  a 
correlation  between  wind  shear  andyH  existed  in  the  upper  levels,  it 
was  necessary  to  examine  the  data  at  times  of  minimum  and  maximum  M 
(i.e.  0>00  OCT  and  1^00  0ST).   Computation  of  the  required  shear  factors 
showed  a  strong  correlation  between  low  shear  values  and  high  M  values 
and  vice  versa  for  all  levels.   Thi3  did  not  controdict  previous  results 
above  600  meters  but  it  was  in  opposition  to  the  relationship  found 
between  shear  and/< below  o00  meters  at  0900  and  2100  GST. 

Since  an  expression  throughout  the  vertical  for  A  in  terms  of 
stability  and  wind  shear  is  a  mo3t  desirable  form,  a  compromise  appeared 
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TABLE  7. 


Stability  factor,  —  — -«l^  ,  versus  height  of  center  of  layer  for 
24  day  average  soundings  at  0300  GOT  and  1^00  GOT, 


0^00  CCT 


La 


,rer 


Height  of 
center  ^ 

^mexers) 


Layer 


Height  of 

center  of 

layer 

(meters) 


1^00  GOT 


Surface  to  950nt 
950nb  to  900ml: 
900mb  to  850m> 
8?0mb  to  SOOmt 
800mb  to  750m 


1 


94 
424 

90s 
1417 
1955 


7498 
,0826 
,0865 
,1065 
,1^1 


Surface  to  950mb 

950mb  to  900mb 

900mb  to  850mb 

850mb  to  800mb 

800nb  to  750mb 


100 

457 
920 

1428 

1965 


,0807 
.19^9 
,1265 

,0969 
,1250 
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to  be  necessary,  i.e.,  a  tv/o  layer  model  with /M  a  function  of  the 
stability  factor  below  600  meters  and/^  a  function  of  the  shear  factor 
above  600  meters.   Now,  from  visual  inspection  of  figs.  1  and  2  at 
1500  GST  and  from  prior  shear  computations  the  shear  can  be  seen  to 
be  very  small  and  practically  constant  throughout  the  layer  surface 
to  1800  meters;  yet  the  value  of  M  (see  fig.  5)  increases  rapidly 
from  600  meters  upwards.   Therefore,  a  second  compromise  from  the 
originally  intended  solution  appeared  to  be  in  order,  i.e.,  above 
600  meters  let  fA  vary  with  the  shear  factor  and  with  height. 
Further,  due  to  the  lack  of  temperature  data  below  600  meters,  it  was 
decided  to  utilize  the  mean^M  (z)  with  respect  to  time  in  an  attempt 
to  more  accurately  describe  the  f4  observed  values  and  also  perhaps  to 
extend  the  coverage  below  o00  meters  to  the  entire  2h   hour  period 
rather  than  merely  the  times  of  the  two  temperature  soundings. 

The  general  approach  used  in  the  solution  of  the  problem  is  as 
follows*. 

1)   Assume  (for  both  layers) 

/i  -  Az)  +  -V* 


(12) 


wherey^-f  is  the  observed  value  of  /M  for  any  point  in  time  and  space, 
yM  (z)  is  the  mean  value  of  all  f^  observed  at  the  level  z,  and  &/4 
is  the  difference  between  the  observed  value  of/^f  and  the  mean  value 
of /M  at  that  height. 

2)   Find  the  best  function  of  height  which  yields  the/M(z) 
distribution. 

5)  Account  for  A/^  ,  the  deviation  from  the  mean, 
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a)  below  600  meters  by  a  function  of  stability  factor, 

b)  above  oOO  meters  by  a  function  of  shear  factor* 

A  function  was  found  in  Section  2  to  represent  the  mean  M  at  any 
height  (see  eq  (7));  however, this  was  based  upon  values  obtained  from 
the  isolines  of/-/  in  fig.  4.   For  this  section  an  attempt  was  being 
made  to  find  the  function  which  yielded  the  computed  values  plotted 
on  fig.  5>  inasmuch  as  they  were  in  part  computed  from  the  vertical 
shear  of  the  horizontal  wind.   A  similar  approach  to  that  used  in 
deriving  eq  (7)  was  in  order  and  yielded  the  following  for /H  (z) 

;q(z)  |  5.5   |z  -4.8  |2'2  +  55  (15) 

where  z  is  height  above  the  ground  in  hundreds  of  meters.   Due  to  the 
relatively  good  fit,  this  equation  was  taken  to  apply  in  both  layers. 

Hex;  A/4  values  were  computed  by  subtracting  the  M  (z)  values 
computed  with  en  (15)  from  the  values  of  M observed  on  fig.  5° 

A  plot  of  A/{  versus  stability  factor  is  shown  as  fig.  12. 
Four  of  the  points  plotted  are  from  observed  data  for  the  period 
(table  7)  and  the  fifth  point  is  an  estimated  value  for  a  height  of 
150  meters  at  1500  GST.   This  value  was  derived  from  12  soundings 
taken  during  the  Great  Plains  Turbulence  Field  Program   5  1  in 
August  and  September  1955  and  under  similar  atmospheric  flow 
conditions.   The  stability  factor  obtained  was  ,0055  n   >  anQ  the 
A/4  corresponding  to  1500  GST  and  150  meters  is*  iHOg  cm"  sec"  . 
The  curve  appearing  on  the  figure  has  been  faired  in  and  does  not 
represent  the  fitted  function  given  below. 

Of  various  functions  of  the  stability  factor  tried,  that 
which  yielded  the  smallest  "least  squares"  error  forAyM  was  the 
hyperbola., 
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-   12.1  (S*  )    4  i  \2O0,5C0  -   74OO(30  )    4   12*      (3@  )^J 


where      (Se  )    =i-  5§    x   107(m'"1)„  (l4) 

Combining  eq  (15)  and  eq  (14)  then  yields  the  desired  form  for 
computing  H  ,  i0e, 

yM  computed  -  5„60  I  z  -  408  2o20  4  186   -  5»6   (3&) 


i    [200,^0  -  74O0(3g>)  4  125     fS^)2]]' 


(15) 


-2 
Here  z  is  in  meters  x  10 

Table  8  shows  a  comparison  of  yM  computed  and  /^  observed  for  the 
layer  below  600  meters „   The  largest  error  occurs  at  2100  GST  and  a 
height  of  424  meters.   It  should  be  noted  (see  figSo  .1,  2)  that  at  this 
point  the  shear  is  relatively  strong  which,  according  to  eq  (3 )  would 
tend  to  yield  a  small  value  f or  /4  <,   A  similar  analysis  applies  for  the 
next  largest  error,  percentwise,  at  0900  C3T  and  4^7  meters „   Here, 
where  the. observed  value  is  greater  than  the  computed  value,  a  low 
value  of  wind  shear  is  observed,,   From  the  preceding  analysis,  it 
appears  that  shear  is  becoming  a  factor  in  the  upper  levels  of  the 
lower  layer o   In  this  transition  zone  from  approximately  400  to  600 
meters,  a  better  representation  of M  might  be  obtained  by  considering 
shear  as  well  as  stability0 

For  the  layer  above  600  meters,  a  graph  was  drawn  showing  height 
versus  shear  factor  and  the  associated  l\M  value  (next  to  each  point) 
needed  to  yield  the  observed  value  of /^  for  the  four  times  shear 
computations  were  made  ,   This  is  shown  as  figo  1J, 
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TABLE  80 
Comparison  of  M     computed  and  M    observed(g  zra"   sec"  ), 

surface  to  oOO  meters. 


Time 

Height  above  Surface 
(meters) 

M   observed 

'        -1   -1 
g  cm  3ec  x 

//computed- 
g  cm"  sec""-'- 

2100 

94 

60 

60 

2100 

h2h 

27 

49 

0900 

1 00 

l4o 

155 

0900 

437 

24 

16 

1500 

150 

4£0 

450 
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As  previously  statei      was  desired  to  account  for  A/M  values  as 
a  function  of  height  an^  shear .   The  approach  to  the  problem  was 

iar  to  that  used  in  section  2„   An  equilateral  hyperbola  was  fitted 
to  the  observed  values  of  shear  and  A/4  at  each  level  (600,  800,  1000 
o<x,ol300  meters)  to  provide  the  A/^  distribution  as  a  function  of 
shear  factor „   The  result  is 


z  - 


ffi 


4-  b 


^N^l 


ll  ..u        "  (16) 


where  l\Mz   is  the  difference  between  the  observed /M  and  M  (z)  at 
height  z  above  the  ground;  c2  and  bz  are  the  constants  determined  by 
the  least  squares  method  at  the  z  level o   Figo  14  shows  c  and  b 
versus  height a   The  curves  appearing  on  the  figure  have  been  faired 
in  and  do  not  represent  the  fitted  function  given  below0 

Next  the  cz  and  bz  values  determined  at  each  height  were  fitted 
with  a  continuous  function  of  height.,   The  functions  of  height 
representing  each  of  these  is 

c(z)  =  ,22    (z)4<  (17) 

b(z)  =  .10   z)^°r  (16) 

-2 
where  z  is  height  above  ground  in  meters  x  10  , 

A  comparison  of  the  cz  and  bz  values  observed  and  the  c(z)  and  b(z) 

values  computed  is  shown  in  table  9«   The  resulting  function  of  shear 

and  height  representing  the  l\m   values  at  600  meters  and  above  is 

.  &      & -(•    )    (^  *  (19) 
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TABLE  9o 
The  constants  c„  and  b_,  versus  c(z)  and  b(z)  at  each  level- 


z 

cz 

c(z) 

bz 

b(z) 

-2 
meters  xlO 

g  cm"  sec"^ 

g  ca  sec"-' 

-1   -1 

< .  cm  sec 

g  cm"  sec 

j> 

314 

509 

-56 

-51 

3 

62J  - 

982 

-119 

-i4o 

10 

56o4 

244o 

-520 

-502 

12 

8596 

4980 

-609 

-574 

14 

7457 

9561 

-870 

-970 

15 

1^,572 

16,115 

-1590 

-1541 

18 

24,147 

25,840 

-2578 

-251s 

41 


and  the  function  representing  M  in  the  upper  layer  is 

/4  =  5.6   (z  -  4„8)2°2    65  +_^22  (z)  4°°  _   -£:)  (z)5 

PW      <», 

A  comparison  of  the  computed  JA  values  and  the  observed  values  is 
shovm  in  Table  10 „ 

The  distribution  of /•'I  computed  indicated  in  Table  10  compares 
fairly  well  with  the  observed  values.   The  salient  features  of  the 
observed^  distribution  are  preserved  by  the  computed  pattern,  i,e,, 
a  maximum  at  each  level  at  time  1500  C3T  (one  exception  at  600  meters) 
and  a  minimum  at  each  level  at  time  0;;00  C3T,   In  most  instances  the 
order  of  magnitude  of  the  observed  and  computed  values  are  the  same 
(five  exceptions  out  of  28  observations).   The  value  observed  at 
0900  GST,  height  o00  meters,  is  approximately  one  fifth  that  computed 
by  shear  and  height  function.   This  discrepancy  may  partially  be 
accounted  for  by  the  higher  stability  factor  observed  in  the  vicinity 
(see  table  7)»   The  negative  values  computed  for  heights  1200  and  l400 
meters  at  0*00  GST  are  primarily  a  result  of  the  curve  fitting  process, 
(Gee  figo  14  and  table  9). 
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TABLE  10 
Comparison  of  computed  coefficient  of  eddy  viscosity ,Mc\g    :m  sec   ) 
and  observed,  UQ    for  four  tines  in  upper  layer  , 


Height 
^meters) 

Mc 

(x0 

A 

(>o) 

A 

H 

A 

<W 

1800 

7607 

(7020) 

455 

(480) 

105 

(157) 

1100 

(528) 

loOO 

6709 

(5900) 

75^ 

(552) 

56 

(97) 

509 

(208) 

lAoo 

5800 

(4680) 

142 

(200) 

-29 

(65) 

177 

(155) 

1200 

1811 

(5080) 

44 

(115) 

-7 

(57) 

120 

(112) 

1000 

856 

(1250) 

41 

(65) 

■7 

(20) 

75 

(81) 

800 

577 

(588) 

24 

W) 

15 

(17) 

46 

(60) 

600 

99 

(142) 

48 

(58) 

26 

(23) 

228 

(42) 

21Z 
150 


05Z 
21c 


osz 

0?C 


0CC 
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ho      Conclusions, 

From  the  approach  followed  in  this  investigation  a  siraple 
representation  of  the  diurnal  variation  of  the  coefficient  of  eddy- 
viscosity  with  height  was  not  apparent. 

It  has  been  shown  that  the  function  of  height  and  time 
representing  the  mean  and  first  two  harmonics  is  inadequate  to  show  the 
salient  features  of  the  observed  distribution. 

In  order  to  represent  the  observed  distribution  of  the  coefficient 
of  eddy  viscosity  by  a  function  of  height  and  time,  it  would  take  a 
relatively  long  series  of  terms,  estimated  to  encompass  at  least  five 
harmonics,  involving  exponential  functions  of  height,  fourth  or  fifth 
degree  polynomial  functions  of  height,  and  trigonometric  functions  of 
angles  dependent  on  time  and  functions  of  height. 

A  two  layer  model  is  developed  in  an  attempt  to  account  for  the 
observed  distribution  of  the  coefficient  of  eddy  viscosity  as  a 
function  of  stability  and  shear  of  the  wind.   It  is  found  useful  in 
the  lower  layer,  and  essential  in  the  upper  layer,  to  include  height 
above  the  ground  as  one  of  the  parameters  determining  the  distribution. 
The  variance  of  the  coefficient  with  stability  in  the  lower  layer 
appears  well  founded,  as  does  that  between  shear  and  the  eddy 
viscosity  coefficient  in  the  upper  layer;  consequently,  a  fair 
representation  of  the  general  characteristics  of  the  observed  diurnal 
variation  of  the  coefficient  of  eddy  viscosity  is  obtained  from  the 
empirical  formulae. representing  the  two  layers. 

Further  investigation  utilizing  more  frequent  and  detailed  wind 

and  temperature  observations  would  be  likely  to  yield  better  results 

« 

and  might  eliminate  the  need  for  the  two  layer  approach. 
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